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ABSTRACT 

We perform a systematic study of hybrid star configurations using several parametrizations of a 
relativistic mean-field hadronic EoS and the NJL model for three-flavor quark matter. For the hadronic 
phase we use the stiff GM1 and TM1 parametrizations, as well as the very stiff NL3 model. In the 
NJL Lagrangian we include scalar, vector and 't Hooft interactions. The vector coupling constant 
g v is treated as a free parameter. We also consider that there is a split between the deconfinement 
and the chiral phase transitions which is controlled by changing the conventional value of the vacuum 
pressure — f2o in the NJL thermodynamic potential by — (fig + SCIq), being SQq a free parameter. We 
find that, as we increase the value of S^to, hybrid stars have a larger maximum mass but are less stable, 
i.e. hybrid configurations are stable within a smaller range of central densities. For large enough SQ Ql 
stable hybrid configurations are not possible at all. The effect of increasing the coupling constant 
g v is very similar. We show that stable hybrid configurations with a maximum mass larger than the 
observed mass of the pulsar PSR J16I4-2230 are possible for a large region of the parameter space of 
g v and SQo provided the hadronic equation of state contains nucleons only. When the baryon octet is 
included in the hadronic phase, only a very small region of the parameter space allows to explain the 
mass of PSR J1614-2230. We compare our results with previous calculations of hybrid stars within 
the NJL model. We show that it is possible to obtain stable hybrid configurations also in the case 
5Qq = that corresponds to the conventional NJL model for which the pressure and density vanish 
at zero temperature and chemical potential. 

Subject headings: stars: neutron — equation of state — PSR J1614-2230 



1. INTRODUCTION 

The recent determination of the ma ss of the pulsar 
PSR J 1614-2230 with 1.97 ± O.O4M bv lDemorest et all 
( 2010), renewed the discussions about the possibility of 
exotic matter being present at the core of neutron stars. 
Since the description of matter at densities beyond nu- 
clear saturation is model dependent, several works have 
explored different aspects of the fact that the maximum 
neutron star mass implied by any equ ation of state must 



exceed t he ma s s of PSR J1614-2230 (Wcisscnbor n et al 



2011afbl. [20121: iBednarek et all 120111: lLastowiecki et al 
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20 1 2t iBonanno fc Sedr akianl 12012^ . Some authors have 
revisited the role of hyperons in the equation of state 
showing that it is possible to construct stiff equations of 
state (EoS) with hyperons that are compatible w ith up- 
to-date hypernuclear data (iBednare k et al.l 120111 ). Oth- 
ers have investigated the role of the vector meson hy- 
peron coupling, going fr om SU(6) quark model to a 
broader SU(3) symmet ry dWeissenborn et al.ll2011bl ) and 
of hyperon potentials (Weis senborn et aLll2012l ) in order 
to determine their impact on the maximum mass of neu- 
tron stars. 

Concerning quark matter, it is known that mod- 
els of strange stars made of absolutely stable 
quark matter satisfy comfortably the new con- 
straint if color-supercondu c tivity is taken into accoun t 
(jLueones fc Horvathl 120031 : iHorvath fc Lugones! 120041 ) . 
However, it is not straightforward to construct mod- 
els of hybrid stars with more than two solar ma sses 
(jBenhar fc CipoUondl201ll: iWeissenborn et al.ll2011al) . A 
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recurrent difficulty is that most hybrid EoS don't have 
at the same time a stable quark matter core and a suffi- 
ciently large maximum mass. For instance, most versions 
of the widely used Nambu-Jona-Lasinio (NJL) model are 
too soft to meet any of the above req uirements. In a re- 
cent work, iBenhar fc Cipollond (|2011l) performed a sys- 
tematic study of the role of the vector and instanton- 
induced terms in the NJL Lagrangian and their connec- 
tion with the properties of hybrid stars. They explored 
a broad region of the parameter space showing that the 
instanton-induced interaction does not affect the stiffness 
of the quark matter EoS, whereas the effect of the repul- 
sive vector interaction is significant. However, according 
to these authors, no values of the corresponding coupling 
constant s allow for the formation of a stable core of quark 
matter (jBenhar fc Cipollond 120111 ) . These conclusions 
are in qualitative agreement with previous results using 
a similar EoS for the hadronic phase but including color 
superconductivity and neglecting the ve ctor interaction 
term in the NJL EoS (jBaldo et al.ll2003D . and with NJL 
models that implement the use of a density-depen dent 
cutoff (jBaldo et al.l[2007t ). Also. lCoelho et all pOlOr ) use 
a SU(2) NJL model with a vector term but they still 
find unstab l e hyb rid stars. Within a different picture, 
iLenzi et al.l (|2010f ) show that it is possible to obtain a 
stable sequence of compact stars with a quark core using 
a NJL model with an ad hoc momentum cutoff that de- 
pends on the baryon chemical potential. However, within 
this approach, the maximum mass i s still below the mass 
of PS R J1614-2230. According to IBenhar fc Cipollond 
(|2011| ). their results are not essentially affected by the 
assumption that the hadronic phase consists of nucle- 
ons only, or by the formation of mixed phases. How- 
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ever, more recent work by Bonanno fc Sedrakianl (pH) 
has succeeded in obtaining very massive stable hybrid 
conhgurations using a three-flavor NJL model that con- 
tains two free parameters: the transition density from 
hadronic matter to quark matter and the vector cou- 
pling of quarks. They show that high-mass stable con- 
figurations with color-superconducting quark cores can 
be constructed if vector interactions are included in the 
quark phase and if a very stiff hadronic equation of state 
is employed (e.g. the NL3 model with hyperons or the 
GM3 model with nucleons). 

In this work we perform an extensive study of hybrid 
star masses using several parametrizations of a relativis- 
tic mean-field hadronic EoS together with a typical three- 
flavor NJL model with scalar, vector and 't Hooft in- 
teractions. Within this approach, the hadronic and the 
quark-gluon degrees of freedom are derived from different 
Lagrangians and the deconfinement transition is associ- 
ated with the point where both models have the same free 
energy. Thus, by construction, chiral symmetry restora- 
tion in the quark model occurs at a chemical potential /i 
that is in general different to the /x of deconfinement. Al- 
though this behaviour is an artefact of the hybrid EoS, 
it might be related with the actual properties of high 
density matter. Present numerical simulations of lattice 
QCD indicate that for small chemical potential /x, the 
deconfinement and chiral transitions coincide, and are 
crossover. At finite ix, a chi ral critical end p oint may 
exist in the phase diagram (|Stephanovl [20051. It has 
been conjectured that, if such a critical end point ex- 
ists, the deconfinement and chiral transitions split from 
one another at that point (|McLerran fc Pisarskil f2007). 
As a consequence, the so called quarkionic phase (con- 
fined but chiral ly symmetric) could be pr esent in the 
phase diagram (jMcLerran fe Pisarskil |2007[ ) . In a sim- 
ilar way, we shall consider that the split between decon- 
finement and chiral symmetry restoration arising in the 
here-employed phenomenological hybrid EoS might rep- 
resent an actual behaviour of matter at low temperatures 
and very high densities. In our model the unknown mag- 
nitude of this hypothetical split will be studied paramet- 
rically by changing the conventional value of the vacuum 
pressure ~Uq which is usually introduced in the grand 
thermodynamic potential of the NJL model in order to 
force a vanishing pressure at zero temperature and /i = 0. 
Instead of — f2o , we shall use a value — (fio + <^o), where 
SQq is a free parameter. Since a change in the value of fio 
has no influence on the fittings of the vacuum values for 
the pion mass, the pion decay constant, the kaon mass, 
the kaon decay constant, and the quark condensates, we 
shall treat it here as a free parameter. Analogously, the 
value of the vector coupling constant g v can be treated as 
a free parameter because the masses of the vector mesons 
are not dictated by chiral symmetry. 

The paper is organized as follows. In Sec. [5] we briefly 
present the here-employed hadronic EoS. In Sec. [3]we de- 
scribe the NJL SU(3) model used to describe the quark 
matter paying particular attention to the role of the pa- 
rameters SD,q and g v . In Sec. |4]we present our results 
and in Sec. [5] we draw our main conclusions. 

2. HADRONIC PHASE 

The relativistic mean-field model is widely used to 
describe hadronic matter in compact stars. In this 



TABLE 1 

Coupling constants for the pa ramet rizations GM1 
IIClendenning He M0SZK0WSKilll991 | ). TM1 ( ISugahara fc TqkiI 
1199411 and NL3 (Ilalazissis et al.II1997I) . M max is the 
maximum mass of a pure hadronic star for matter 
composed by nucleons and electrons. 



Set 


GM1 


TM1 


NL3 


m a (MeV) 


512 


511.198 


508.194 


m„ (MeV) 


783 


783 


782.501 


m p (MeV) 


770 


770 


763 




8.91 


10.029 


10.217 




10.61 


12.614 


12.868 


9p 


8.196 


9.264 


8.948 


b 


0.002947 


-0.001506 


0.002055 


c 


-0.001070 


0.000061 


-0.002651 


Mmax 


2.32 


2.18 


2.73 



paper we adopt the fol lo wing standard Lagrangian 
jlBoguta fc Bodmerl 119771 : IGlendenning fc Moszkowskil 
119911): 



B^ 



:9 P BT.p^) 



(m B - 9aB(r)]ipB + - {d^ad^a - m a a ) 



(1) 



for matter composed by (i) nucleons and electrons, and 
(ii) the baryon octet and electrons. Leptons L are 
treated as non-interacting and baryons B are coupled 
to the scalar meson <r, the isoscalar- vector meson w M 
and the isovector- vector meson p p . For more details 
about the EoS obtained from the above Lagrangian the 
reader is referred to e.g. iLugones et al.l (|2010( ) and ref- 
erences therein. There are five constants in the model 
that are fitted to the bulk propertie s of nuclear mat- 
ter ([Glcndcnni ng fc Mo szkowski 199l|). In this work we 
use three different parametrizations shown in Table [TJ 
For all parametrizations we use a composition of nucle- 
ons and electrons. For the NL3 parametrization we con- 
sider also the case with the baryon octet and electrons. 
The parametrization for the hyperon coupling constants 
is gujA/g^N = g^v/g^N = 0.6666, g u s/guN = 0.3333, 

gak/gaN = 0.6106, q„y Jq„N = 0.4046, Q^./q^ N = 



0.3195 and g pi /g P N = 1 (|Chiapparini et alfl2009l ). At 



low densities we use the Baym, Pethick and Sutherland 
(BPS) model (jBavm et al.lll971l ). 

3. QUARK PHASE 

3.1. The model 

To describe the quark matter phase we use the SU(3) 
NJL model with scalar-pseudoscalar, isoscalar-vector 
and 't Hooft six fermion interaction. The Lagrangian 
density of the model is: 
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+ fft {det#(l + 75)^] + det$(l - 75)$}, (2) 

where i/j — (u, d, s) denotes the quark fields, X a (0 < a < 
8) are the U(3) flavour matrices, rh = diag(TO„, mj, m s ) 
is the quark current mass, and g s , g v and gt are coupling 
constants. 

Notice that we have not included a diq uark interac- 
tion t erm in the Lagrangian. As shown bv lRiister et al.l 
(2005), color superconducting phases (at T = 0) are 
favoured in the regime of strong diquark coupling, 
9d/9s ~ 1- However, in the regime of intermedi- 
ate diquark coupling strength, g c i/gs — 3/4, color 
superconductivity appears only above a chemical po- 
tenti al » ~ 3 x 440 M eV = 1320 MeV (see Fig. 
1 of iRiister et alj (|2005|) ). For weaker diquark cou- 
pling, color superconductivity is shifted to very large 
densities that are not present at neutron star cores. 
Since the case of st rong diquark coupling has al ready 
been considered bv iBonanno fc Sedrakianl (|2012l ) and 
iPagliara fc Schaffner-Bielichl (12008) , we shall focus here 
in a case where color superconductivity is negligible. 

The mean-field thermodynamic potential density for 
a given baryon chemical potential \x at T = 0, is given 

by 



A o. 



p dp 
~2n r 



p2 + Mf + 2g s J2{Ml 



vNcJ2 [ 

-2g v J2(^^+4gt(uu)(dd)(ss) 

i 

mX- f kF 'P 2 dp 



(3) 



where the sum is over the quark flavor {% = u,d,s), 
the constants 77 = 2 and N c = 3 are the spin and 
color degeneracies, and A is a regularization ultravi- 
olet cutoff to avoid divergences in the medium inte- 
grals. The Fermi moment of the particle i is given by 
k Fl = 9(/j,* - Mi)yJ~(pL* 2 - Mf), where pi* is the quark 
chemical potential modified by the vectorial interaction, 

In this work we consider the follow ing set of parameters 
(jKunihirol fl989t IR.uivo et al.l fl999h : A = 631.4 MeV , 
g s A 2 = 1.829, g t A 5 = -9.4, m u = m d = 5.6 MeV, m s = 
135.6 MeV, in order to fit the vacuum values for the 
pion mass, the pion decay constant, the kaon mass, the 
kaon decay constant, and the quark condensates: — 
139.0 MeV, = 93.0 MeV, m K = 495.7 MeV, f K = 
98.9 MeV, (to) 1 / 3 = (ui) 1 / 3 = -246.7 MeV, (ss) 1 / 3 = 
—266.9 MeV. The value of the vector coupling constant 
g v is treated a free parameter because the masses of the 
vector mesons are not dictated by chiral symmetry. In 
order to obtain the equation of state, we assume that 
matter is charge neutral and in equilibrium under weak 
interactions. 

3.2. fig as a free parameter 
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Fig. 1. — (a) Pressure as a function of the chemical potential for 
different values of the parameter SQq. (b) Pressure of the decon- 
finement phase transition as a function of 8Q0 for different values 
of the coupling constant g v . Notice that a small change in <5f2o can 
produce a significant change in the pressure of the phase transition. 

The conventional procedure for fixing the f2o term 
in Eq. ^ is to assume that the grand thermody- 
namic potential 51 must vanish at zero /j, and T. For 
the above quoted parametrization, this assumption leads 
to the value fl = 5076.2 MeV fm -3 . Nevertheless, 
this prescription is no more than an arbitrary way to 
uniquely determine the Eo S of the NJL model without 
any further assumptions (jSchertler et al.l Il999h . Fur- 
thermore, in the MIT bag model for instance, the 
pres sure in the vacuum is non-vanishin g. In view of 
this, IPagliara fc Schaffner-Bielichl (|2008l) adopt a differ- 
ent strategy. They fix a bag constant for the hadron- 
quark dcconfinement to occur at the same chemical po- 
tential as the chiral phase transition. This method leads 
to a significant change in the EoS with respect to the 
conventional procedure. 

The connection between the chiral and the deconfine- 
ment phase transitions along the QCD phase diagram 
has received considerab l e atte ntion in recent years (see 
iFukushima fc Hatsudal (| 2 1 ID and references therein). 
For zero chemical potential, lattice results sh ow that 
both t ransitions occur at the same temperature ([Karschl 
120021: lLaermann fc Philipsenl I2003D . At finite baryon 
chemical potential, this coinc idence is an open question 
(|Fukushima fc Hatsuda 1201 ID . However, since a chiral 
critical end poin t may exist in the plane of T and p 
(|Stephanovl I2005D . it has been conjectured that the de- 
confinement a nd chiral transitions split fro m one another 
at that point (Mc Lerran fc Pisars 12001 . As a conse- 
quence, a confined but chiral symmetric phase, which 
is called quarkyonic phase can exist in the high baryon 
density region. Since this conjecture is based on argu- 
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Fig. 2. — (a) Pressure as a function of the baryon number density 
in units of the nuclear saturation density po (we assumed po = 
0.17 fm — 3 ). (b) Mass of hybrid stars as a function of the central 
mass-energy density e c . We use <5f2rj = and different values of g v . 

ments that are valid in the large N c limit, it is not clear 
whether this quarkyonic phase can exist in the real QCD 
phase diagram. In the present paper we are not mod- 
elling the quarkionic matter because our confined phase 
is described by a hadronic model and the chiral transi- 
tion is restricted to the quark phase. However, we may 
explore the above possibility of having chiral restoration 
and deconfinement occurring at different densities. To 
this end, we shall substitute ilo in Eq. (j3j by the new 
value Qq + 5£Iq, where <5f2o is a free parameter: 

n — ■+ O + 8n in Eq. ©. (4) 

With this change, the thermodynamic potential f2 can 
be non- vanishing at zero fj, and T, and the fi of the de- 
confinement transition can be tuned. Clearly, <5f2o nas 
a minimum value because the phase transition cannot 
be shifted to a pressure regime where the NJL model 
describes the vacuum. That is, we fix a minimum 
limit to SQq for which the phase transition occurs at 
the chiral symmetry restora t ion p oint as performed by 
iPagliara fc Schaffner-Bielichl (|2008t l. On the other hand, 
there is no maximum value in principle for S^Iq, since 
the phase transition can be shifted to arbitrarily large 
pressures. 

In order to illustrate the dependence of the EoS on the 
new parameter dflo we depict in Fig. Q] the pressure as 
a function of the chemical potential for different values 
of SQq and the pressure of the deconfinement transition 
Ppht as a function of <5f2rj. Notice that a small change in 
the value of <5fio may result in a significant modification 
of the phase transition density, and consequently, in a 
very different hybrid EoS. 

4. RESULTS 

We have solved the Tolman-Oppenheimer-Volkoff 
equations for spherically symmetric and static stars in 
order to investigate the influence of g v and <5f2o on the 
maximum mass of hybrid stars. 

In Figs. [2] and [3] we show the EoS for some specific 
parametrizations and the corresponding stellar configu- 
rations in a diagram of mass M versus central energy 
density e c . The plateaus represent the hadron-quark 
phase transition as a consequence of a first order Maxwell 
construction. In the left panel of Fig. [2] we note that as 




Fig. 3. — Same as Fig. [2]but adopting g v /ga = 0.2 and different 
values of 5f2o (labels for <5f!o are in MeV fm" 3 ). 

we increase the value of the vector coupling constant g v 
so does the density of the phase transition, and there- 
fore, the hybrid star has a smaller quark core and a 
larger hadronic contribution. This leads to larger maxi- 
mum masses because the hadronic EoS is stiffcr than the 
quark EoS. At the same time, there is a larger density 
jump between the two phases, which tends to destabi- 
lize the configuration. Due to these two effects, together 
with the fact that the vector term stiffens the NJL EoS, 
models with a larger g v give larger maximum masses but 
have stable quark cores within a smaller range of central 
densities (see right panel of Fig. [5]). In the left panel of 
Fig. [3] we show the effect of changing the magnitude of 
the shift between the deconfinement and the chiral phase 
transitions. As we increase SQq from negative to positive 
values, we increase the density of the phase transition as 
well as the density jump between the two phases. How- 
ever, the NJL EoS becomes slightly softer because there 
is a larger contribution to the EoS of the regime with a 
partially restored chiral symmetry. Since the latter effect 
is not so strong, the impact of increasing 5Qq is analo- 
gous to increase g v , i.e. models with a larger SQq result 
in larger maximum masses but the quark cores are sta- 
ble within a smaller range of central densities (see right 
panel of Fig. [3]). 

In Fig. |4l background colors represent the maximum 
mass of hybrid stars for different values of g v and SQq. 
We used the hadronic EoS of Table [T] with nucleons and 
electrons. Within each panel we show contour lines indi- 
cating specific values of the maximum mass. The black 
line represents the limit between parameters that allow 
for stable hybrid stars and those that always give un- 
stable hybrid stars. The value 1.97M Q c orresponding to 
the observed mass of PSR J1614-2230 (jDemorest et all 
2010; is shown with a red dashed line. An interesting 
feature of Fig. QJ is that large masses are situated on the 
right-upper corner but stable configurations are located 
on the left-lower corner of the figure (or left side of the 
figure in the case of NL3). This clearly illustrates the 
difficulty of obtaining stable hybrid stars with arbitrar- 
ily large masses. Concerning the effect of the hadronic 
model we see that stable hybrid stars have higher values 
of the maximum mass for the stiffer hadronic EoS. 

The observed mass of PSR J1614-2230 can be ex- 
plained by parameters within the large region located 
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Fig. 4. — Background colors represent the maximum mass of 
hybrid stars for different parametrizations of the NJL model (i.e. 
different values of g v and <5f2o). fn each panel we use a differ- 
ent hadronic EoS (without hyperons): (a) GM1, (b) TM1 and 
(c) NL3 (see Table [TJ- Notice that the color scale is different 
for each panel. The solid contour lines indicate specific values of 
the maximum mass. The black solid line represents the boundary 
between parametrizations that allow for stable hybrid stars and 
parametrizations that do not. The red dashed line indicates the 
value 1,97-Mp) correspondin g to the observed mass of PSR J1614- 
2230 (Dcmorcst ct al. 2010). The region between the red dashed 
line and the solid black line allows to explain the mass of PSR 
J1614-2230. 

between the red dashed line and the solid black line in 
each panel of Fig. SJ However, a hypothetical future ob- 
servation of a neutron star with a mass a ~ 10% larger 
than the mass of PSR J1614-2230 will be hard to explain 
within hybrid star models using the GM1 and TM1 EoS 
(see panels (a) and (b) of Fig. [4| and will require a very 



,60 

~~~> 

60 




-16 -14 -12 



8Q Q [MeVfm 



Fig. 5. — Same as panel (c) of Fig. 4 but for the NL3 model with 
hyperons. Hybrid stars are not possible for the set of parameters 
within the white region. Only a very small region near the upper- 
left corner of the coloured region allows to explain the mass of PSR 
J1614-2230. 

stiff hadronic model such as NL3. 

The effect of hyperons is shown in Fig. [5] where we 
consider the NL3 parametrization with the inclusion of 
the baryon octet. Compared with the case without hy- 
perons, the maximum mass values are altered by a few 
percent. This follows from the fact that the deconfine- 
ment phase transition occurs at relatively low densities, 
i.e. in regions where the baryon octet has a minor contri- 
bution. Nevertheless, hyperons have a large effect in the 
possibility of finding stable hybrid stars of large enough 
mass. When we increase the value of g v or increase SQo, 
the deconfincment transition is shifted to larger densi- 
ties and the hadronic EoS with hyperons tends to be 
favoured almost everywhere in the star. Above a cer- 
tain limit there is no deconfinement transition at all, i.e. 
the star is always hadronic (see white region in Fig. 5). 
As a consequence, the mass of PSR Jf6f4-2230 can be 
attained for models within a very small region of the pa- 
rameter space located near the upper-left corner of the 
coloured region in Fig. 5. 

5. DISCUSSION AND CONCLUSIONS 

In this work we performed a systematic study of hy- 
brid star configurations using a relativistic mean-field 
hadronic EoS and the NJL model for three-flavor quark 
matter. For the hadronic phase we used the stiff GM1 
and TM1 parametrizations, as well as the very stiff NL3 
model. In the NJL Lagrangian we included scalar, vector 
and 't Hooft interactions. The vector coupling constant 
g v was treated as a free parameter. We also considered 
that there is an arbitrary split between the deconfinc- 
ment and the chiral phase transitions. This split can be 
adjusted by a redefinition of the constant parameter fio 
in the NJL thermodynamic potential; i.e. by making the 
replacement fl^ — > fig + Sflo in Eq. ([3]), where SCIq a 
free parameter. We find that, as we increase the value of 
5£Iq, hybrid stars have a larger maximum mass but are 
less stable (i.e. hybrid configurations are stable within a 
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smaller range of central densities). For large enough SQq, 
stable hybrid configurations are not possible at all (see 
Fig. 3). The effect of increasing the coupling constant g v 
is very similar (see Fig. 2). These effects are clear in Fig. 
IU where we show the maximum mass of static spherically 
symmetric stars in the parameter space of g v and SQq. 
The larger masses are situated on the right-upper cor- 
ner of the diagram, where both g v and SQq are larger. 
On the other hand, stable configurations are placed on 
the opposite part of the diagram; i.e. on the left-lower 
corner for the GM1 and TM1 models and on the left 
side for the NL3 model. As a consequence, stable con- 
figurations with a maximum mass compatible with PSR 
J1614-2230 are located halfway, specifically, between the 
red dashed line and the solid black line of Fig. 4. The 
effect of the hadronic model (with nucleons only) is also 
clear from Fig. 4, where we see that stable hybrid stars 
have higher values of the maximum mass for the stiffcr 
hadronic EoS. The main effect of hyperons is that they 
preclude the deconfincment transition in the region of 
the parameter space that allows large maximum masses 
(see Fig. 5). Only a very small area near the upper- left 
corner of the coloured region of Fig. 5 allows to explain 
the mass of PSR J1614-2230. 

It is worth summarizing the main assumptions and 
results of recent wor k using the NJL mod el to de- 
scribe hybrid stars. Schertlcr et al.l (|l999) use sev- 
eral parametrizations of an extended relativistic mean 
field model to describe the hadronic phase with hyper- 
ons. For the deconfined quark phase they use the NJL 
model with three flavors, including scala r and 't Hooft 
term s, and using the RHK parameter set (Rchbcrg ct al. 
1996). Their conclusion is that typical neutron stars with 
masses around 1.4 solar masses do not possess any de- 
confined qu a rk ma tter in their center. More recently, 
iBaldo et all (|2007| ) use a SU(2) NJL model with the 
RHK parametrization and a cut-off that depends on 
the chemical potential. For the hadronic phase they 
adopt a nucleonic equation of state obtained within 
the Brueckner-Bethe-Goldstone approach using the Ar- 
gonne v\g two-body potential, supplemented by the Ur- 
bana phenomenological three-body force. They are not 
able to obtain stable hybrid stars. iCoelho et al.l (|2010f ) 
use a SU(2) NJL model with the RHK parametriza- 
tion but adding a vector term together with the GM1 
parame trization of a relativistic mean field hadronic 
model (jGlendenning fc Moszkowskilll991l ). The density 
of the phase transition incr eases with q v but th e y stil l 
find unstable hybrid stars. iBenhar fc Cipollo"nel (|2011[ ) 
work with SU(3) and the RHK parametrization. They 
use several values of g v and they also vary the coupling 
constant of the 't Hooft term. For the hadronic phase 
they use a phenomenological Hamiltonian including the 
Argonne v' e nucleon-nucleon potential. They conclude 
that no values of the corresponding coupling constants 
allow for the formation of a stable core of quark matter. 

All the above versions of the NJL model use the con- 
ventional procedure of imposing that the pressure and 
density must vanish at zero temperature and chemical 
potential. However, a differ e nt pr escription is used by 
IPagliara fc Schaffner-Bielichl (|2008l ). First, they intro- 
duce a hadronic EoS and compute the transition to quark 
matter by a Maxwell construction. To fix the bag con- 
stant they assume that deconfinement occurs at the same 



chemical potential as the chiral phase transition, i.e. they 
require that the pressure of quark matter is equal to the 
pressure of the hadronic matter at the critical chemi- 
cal potential for which chiral symmetry is restored. The 
bag value obtained with this assumption is the lowest 
possible value for the bag constant in the NJL model 
because it allows to use the NJL EoS just starting from 
the chemical potenti al of the chiral phase transition. Us - 
ing this procedure, IPagliara fc S chaffner -Bielichl (2008) 
computed the equation of state of quark matter within 
the NJL model by including effects from the chiral con- 
densates, the diquark coupling pattern, and a repulsion 
vector term. They find that hybrid stars containing a 
CFL core are stable but the maximum mass is ~ 1.8M©, 
i.e. incompatible with PSR J1614-2230. More recently, 
IBonanno~fc Scdrakian (2012) use a NJL model supple- 
mented by the 't Hooft and vector interactions and con- 
sider the 2SC and CFL color superconducting phases. 
For the hadronic phase they use a relativistic mean field 
model and adopt the NL3 parameterization with hyper- 
ons and the GM3 parameterization with nucleons only. 
In both cases they are able to obtain maximum masses 
above the mass of PSR J1614-2230 because they use the 
non-conventional procedure of treating the density of the 
deconfinement transition as a free parameter. 

Our analysis is related to that of lBonanno fc Scdrakian 
(2012) because we can control the density of the phase 
transition via the parameter SQ,q in the thermody- 
namic potential. However, since our parameter space 
is constructed in terms of Sttg and not in terms of 
the de confinement density as in iBonanno fc Scdrakian 
(|2012f) . the connection with conventional NJL mod- 
els is more transparent in our case. Additionally, our 
procedure includes the non-conven tional prescription of 
IPagliara fc Sch affner-Biclicbl (|2008l ) as a special case. It 
is also worth noticing th at we use t he pa r ametrization 
of the NJL EoS given by iKunihirol (|1989t ); iRuivo et al.l 
(1999) while the above a uthors us e the R KH one which 
is somewhat softer (see Buballa (2004) for more de- 
tails on the parametrizations). Another difference is 
that they consider quark matter in the 2SC and CFL 
phases, whereas we don't consider color superconduc- 
tivity in the quark phase. In this sense, these works 
are complimentary because IPagliara fc Schaff ner-Bielichl 
( 2008) ; IBonanno fc Scdrakian (2012) work in the regime 
of strong diquark coupling {ga/gs ~ 1) where color su- 
perconductivity is strongly favoured, whereas we work 
in the weak diquark coupling regime (gd/gs — > 0) for 
which color superconductivity is shifted to very large 
densities tha t are n ot present at neutron star cores (see 
iRiister et all pOOl for more details). 

In summary, our results show that hybrid configura- 
tions with maximum masses equal or larger than the ob- 
served mass of PSR J1614-2230 are possible for a signif- 
icant region of the parameter space of g v and SQq pro- 
vided a stiff enough hadronic EoS without hyperons is 
used. It is also worth highlighting the fact that we can 
obtain compact stars with stable quark cores without 
having to perform any modification to the NJL model 
(i.e. setting 6 flp = 0). Thi s is in c on trast with the results 
obtain ed b y IBaldo et all (120071) . IBenhar fc Cipollond 
(|2011D and ICoelho et al.l (|2010D who are unable to ob- 
tain stable hybrid star configurations with the NJL 
model. The difference arises because these authors 
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use the softer parametrization of iRehberg et all (119960 
while we use the stiffer para metrization of iKunihirol 
(|1989D and iRuivo et all (|1999D : i.e. the use of the lat- 
ter parametrization allows to reproduce the mass of the 
pulsar PSR J1614-2230 without using exceptionally stiff 
parametrizations of the hadronic EoS and keeping the 
conventional procedure for fixing firj. It is also interest- 
ing to note that we have typically Sflo/Qo ~ 0.1%, i.e. 
small departures of the conventional £Iq are sufficient to 
obtain massive enough stable hybrid stars. Finally, we 
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